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We give a new and powerful method for the direct solution of circuit 
design problems. The method begins with a prespecified topology and some 
or all elements undetermined in value. The designer imposes on the circuit 
any desired set of node-pair voltages, branch currents, or driving point 
and transfer immittances. Values of circuit elements that satisfy the con- 
straints are directly calculated. This direct method of solution avoids the 
usual iterative analysis-optimization schemes, reducing computer times by 
up to three orders of magnitude. 

A linear set of design equations is formulated by choosing undetermined 
element currents and node voltages as the variables. Singular elements are 
introduced to impose the desired constraints. Inequality as well as equality 
constraints are permitted. Element values are determined from the solution 
of these equations. In this paper we emphasize our method of solution in 
relation to dc networks. 

I. INTRODUCTION 

The most significant advances made in computer-aided circuit de- 
sign have been in analysis programs. The designer can now choose 
from among several general purpose programs that program which 
most nearly suits his particular needs. In designing a circuit to meet 
a given set of requirements, the usual approach has been to use 
analysis programs in some optimization scheme. Through an iterative 
process, carried out by the machine, the man, or a man-machine inter- 
action, a final design is reached. The approach presented here pro- 
vides a direct solution, and does not rely on such iterative schemes. 

The method is most fertile in the area of active network design, 
where one often wishes to choose element values in a specified topology 
in order to meet some set of requirements. The method has been 
applied to a number of design problems of current interest including 
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biasing direct coupled transistor circuits; designing transistor ampli- 
fiers for specified midband gain, input, and output impedances; and 
simultaneously realizing several specified impedance or admittance 
parameters of a network. 

In the design of electronic circuitry, one usually wishes to imbed 
passive elements into a network containing active devices, and to 
determine the required passive element values. Therefore, this paper 
deals with the determination of element values in a prespecified 
topology for which a given performance is required. Two new ele- 
ments, a voltage forcing element (VFE) and current forcing element 
(CFE), are introduced in order to constrain network voltages and 
currents. These elements may be realized with independent voltage 
and current sources, and the nullator, a somewhat "pathological" 
element used in theoretical network studies. 

The method of singular imbedding places the VFE's and CFE's 
in a network to constrain the desired variables. The terminal voltage- 
current behavior of the variable elements is not specified. Instead, 
the constraints imposed upon the network by the VFE's and CFE's 
are used to determine allowed voltage-current relations for the variable 
elements. The formulation remains linear in these variables. The last 
step involves determining the element values through Ohm's law 
once the allowed voltage-current relations are known. 

By appending the original set of equations with a set of inequality 
constraints, it is possible to restrict the range of element values in 
the solution. For example, realizations employing only element values 
between specified lower and upper bounds are possible. For simplicity, 
only the case of linear dc networks are illustrated. Extensions of 
the method to ac and nonlinear design are considered elsewhere. 

II. A NEW APPROACH 

To understand the philosophy of this new approach to design, con- 
sider the train of events in realizing a set of requirements with elec- 
tronic circuitry. Since the choice of topology is better handled by the 
man than the computer, we will assume some specified topology in 
which some or all of the element values are to be chosen to meet the 
given criteria. For example, in designing transistor circuitry it is nec- 
essary to choose some resistance values to properly bias the transis- 
tors. Similarly, one must often choose element values to give a desired 
voltage gain, driving point impedance, transfer impedance, or similar 
network function. 
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The invariant feature in all of these problems is that a set of net- 
work currents and voltages, or their ratios, has been constrained. The 
design problem is to find any set of element values consistent with 
these constraints. If the problem is posed with sufficient freedom, 
many sets of element values may exist consistent with the imposed 
constraints. Conversely, if the problem is posed with insufficient free- 
dom, inconsistent equations arise and there is no solution. 

If one can find a general method of imposing these network con- 
straints, and can simultaneously monitor the voltage-current relations 
these constraints force at the terminals of the variable elements, then 
indeed a direct solution to many computer-aided design problems will 
have been found. 

Before proceeding, however, consider a very simple example of 
how one might presently handle the design problem and the diffi- 
culties that would ensue. Suppose in the network of Fig. 1, one wishes 
to choose G\ and G 2 such that V is constrained to be 0.1 volt. A set 
of nodal equations may be written : 
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The first step involves a transformation of coordinates 
the desired quantities appear explicitly in the equations. In 
this will necessitate using hybrid parameters. For this case, 
lowing transformation might be used : 
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Fig. 1 — Simple design problem. 
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and the current-law equations become 
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Substituting the constant V = V\ = 0.1, the set of equations becomes 
0.1(1 + <?,) + Vl = 1 
-0.1(7, + (? 2 F 2 f = 0. 

Thus, even if one is successful in finding a transformation to a 
basis that includes the variables that are constrained, the result is 
usually a set of nonlinear equations in the network elements and 
voltage variables. Solving this set of nonlinear equations for the 
unknown voltages and element values is extremely difficult. A method 
of handling this difficulty has been suggested, involving the use of 
optimizing techniques to vary element values until the network vari- 
ables take on their desired values — in this case V = 0.1 volt. 1 While 
this is a useful approach, it has several disadvantages. First, it is 
time consuming since many iterations are required for convergence. 
Second, local minima, or lack of sufficient numerical accuracy, may 
prevent convergence to a correct solution. Finally, although an infinity 
of sets {G x , G 2 ) exist to satisfy the given constraints, the optimization 
yields only one of these sets. 

With these difficulties in mind, let us repeat the philosophy of 
design presented here. We first determine how the requirements con- 
strain network currents and voltages. We then force these currents 
and voltages to take on the desired values. Finally, we determine the 
effect of such constraints upon the voltage-current relations at the 
terminals of variable elements. These v — i relations then determine 
the values of the variable elements. 

III. NETWORK CONSTRAINTS 

The common feature of all network synthesis problems is that they 
require some specified relation between some voltages and currents in 
the network. For example, synthesis of a given driving point im- 
pedance constrains the ratio of a port voltage to the current at that 
port. Synthesis of a transfer impedance constrains the ratio of a port 
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voltage to the current at a different port. A specified voltage or cur- 
rent gain constrains the ratio of two port voltages or port currents, 
respectively. Indeed the synthesis of entire network matrices is a 
combination of such constraints. Similarly, the static design problem 
in electronic circuits involves fixing certain branch currents and 
branch voltages. For example, one usually wishes to bias a transistor 
for a given collector current and collector-emitter voltage. Resistance 
values are chosen consistent with these constraints. 

It is essential to demonstrate a method for constraining voltages 
and currents in a network. The required constraints are shown in 
Fig. 2. We introduce two new elements, a current forcing element, 
CFE(7 ), and a voltage forcing element, VFE(V ), which will be real- 
ized with more conventional elements shortly. We want the CFE(J ) 
to be such that it constrains the current through branch ; to be I , 
without otherwise affecting the behavior of the network. We want the 
VFE(V ) to be such that it constrains the voltage across branch ; 
to be V without otherwise affecting the behavior of the network. 

In discussing the properties of the CFE and VFE, we use the concept 
of admissible or allowed pairs of voltage and current variables (v, i). 2 
The set of voltage-current pairs that a system N allows can be used to 
completely describe that system. 3 For example, let the system under 
consideration, N R , consist of a single resistor of value R. Then the 
system is completely described by its allowed terminal voltage and cur- 
rent pairs; namely, (Ri, i) e N R . Similarly, a capacitance of value C, 
denoted N c , is completely described by its allowed pairs (v, d(Cv)/dt) t 
N c . 

We now define the CFE(/ ) and VFE(7 ) in terms of their allowed 
pairs. 

Current forcing element (/„) : 

(0, /„) e N CFEUo) . (7) 

Here we postulate an element which allows no voltage drop across 
its terminals, and passes only a specified current I . 

Next, we postulate an element which allows only a fixed voltage 
V to exist at its terminals, and passes no current. 

Voltage forcing element (V ) : 

(7 o ,0) zN VFElVo) . (8) 

Figure 2 makes clear the use of these elements in constraining net- 
work variables. In Fig. 2a, the current in branch ; is forced to be 
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Fig. 2 — Network constraints, (a) Branch J current constrained by current 
forcing element (CFE) ; (b) Branch J voltage constrained by voltage forcing 
element (VFE). 

I by inserting a CFE in series. Since the CFE(7 ) allows no voltage 
to exist across its terminals, its presence affects Kirchhoff's current 
and voltage laws only to the extent that branch ; current is con- 
strained to be I n . Notice that this would not be the case had we in- 
serted a current source in series with branch ;'. The current source 
would allow some voltage to exist between its terminals which would 
have been included in Kirchhoff's voltage law equations. Thus, a 
current source of value I would not only constrain branch ; current 
to be I , but would also introduce a new degree of freedom, namely, 
the voltage across the current source. 

Similar reasoning can be applied to Fig. 2b. Here a VFE(7 ) is 
applied across branch j to constrain that voltage to be V . Since the 
VFE(7 ) passes no current, Kirchhoff's laws are affected only to the 
extent that branch ; voltage is now constrained to be V . The net- 
work cannot respond with a new degree of freedom, as it could if a 
voltage source were placed across branch ; and thus allowed to in- 
troduce a new current variable in Kirchhoff's current law equations. 
It should be noted that the VFE(7 ), can be placed between any 
two nodes to constrain the voltage between those nodes; it need not 
be placed across a branch. 

By using current sources and voltage sources in conjunction with 
VFE's and CFE's, current- voltage ratios may be constrained. For 
example, in Fig. 3a, 

Zi - I± - %L 
h I I 

In Fig. 3b, 
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Fig. 3 — Methods of constraining current-voltage ratios, (a) Impedance forcing 
element [IFE(Z)]; (b) admittance forcing element [AFE(Y)]. 

Thus we are constraining the network N to have, in the first case, 
a driving point impedance Z, and in the second case, a driving point 
admittance Y. The configurations used to constrain impedances or 
admittances will be denoted impedance forcing elements, IFE(Z), 
and admittance forcing elements, AFE(Y). Notice that IFE's and 
AFE's are composed of CFE's, VFE's, and independent sources. They 
are useful in constraining a network to have a desired driving point 
impedance or admittance. 

We already mentioned that VFE's and CFE's could be realized in 
terms of existing elements. The necessary elements are the ideal cur- 
rent source, the ideal voltage source, and the nullator, a somewhat 
"pathological" network clement introduced by Tel logon. * Returning to 
the allowed pair concept, the nullator is defined to be a two-terminal 
element for which the only allowed voltage-current pair is (0, 0). It 
can be looked upon as a simultaneous open and short circuit, since it 
allows only zero voltage at its terminals and passes no current. 

From its definition, one could not hope to physically realize and 
isolate such a device. However it's characteristics may be observed at 
the input to an operational amplifier imbedded in a feedback net- 
work, where the input is at a virtual ground (short circuit) and yet 
passes no current (open circuit). The nullator is represented sche- 
matically in Fig. 4. 

By appropriate connections of voltage sources, current sources, 
and nullators, the VFE's and CFE's may be realized as in Fig. 5. 
Remembering that the nullator passes zero current and has zero volt- 
age across its terminals, the equivalents of Fig. 5 becomes clear. In 



Fig. 4 — Schematic representation of nullator. 
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Fig. 5 — Equivalent circuits for VFE and CFE using nullators. 

Fig. 5a, the terminal voltage must be V , and since no current exists 
in the element the combination voltage source and nullator is by 
definition a VFE(V ). In Fig. 5b, a current I exists at the terminals 
but no voltage drop exists across the terminals. Thus by definition, 
the combination current source and nullator is a CFE(7 ). 



IV. ADDING FREEDOM TO THE NETWORK 

In the previous section, we placed constraints on the network that 
would generally lead to a set of inconsistent equations if all the 
elements were also specified. However, if some network elements are 
variable, we can determine how the constraints affect the voltage- 
current relations at the variable element terminals, and then choose 
variable elements in such a way as to be consistent with these v — 
i relations. 

We propose two methods of characterizing the variable elements. 
First, since the element is variable, we can ascribe no functional rela- 
tion between the voltage and current of that branch. This is handled 
in writing the nodal equations for the network by explicitly adding 
the currents through variable elements into the equations, rather than 
first transforming them into voltage variables through a functional 
relation of the form 

it = Y b v b (11) 

where the b implies the variable refers to some branch. The nodal 
equations are of the form 



[Y f ]V] = I,] + [C]I], 



(12) 



where 



V] is an n-vector of node voltages. 

I s ] is an n-vector of forcing currents at each node. 
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[Y f ] is the n X n nodal admittance matrix of the fixed portion of the 

network. 
I] is an r-veetor of unknown currents through variable elements (r 

is the number of variable elements) . 
[C] is the n X r node cutset matrix for the graph of variable elements. 
I] and V] are both vectors of network variables, and may be combined 
by matrix partitioning as 

"" I 



[-C | Y f ] 



V 



= Is]. 



(13) 



Equation (13) describes a network in which some element values 
can be chosen to meet the given constraints. In the remainder of this 
paper, we combine the added degrees of freedom given by the variable 
elements in equation (13) with the constraints imposed by the CFE's 
and VFE's. All networks, satisfying the VFE and CFE constraints 
and the specified topology, with be generated. 

A simple example will help clarify these concepts. Figure 6 is the 
network of Fig. 1, with the 1-ohm resistor replaced by a known resis- 
tance of R ohms. Currents /x and I* are those carried by the variable 
conductances G x and G 2 , respectively. The set of nodal equations is 
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Rearranging into the form of equation (13) , 

h 



(15) 



1 i 1 

i 

L-i i ! o oj v, 

From this example, the method of generating equation (13) should 
become clear. 




Fig. 6 — Simple design problem. 
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The second approach useful in dealing with variable elements in a 
network is the introduction of another pathological element, the nora- 
tor, shown in Fig. 7, also introduced by Tellegen.* The norator is a 
two-terminal element with allowed pairs (v, i), with v and i independ- 
ent and arbitrary. Thus, any voltage and current may appear across 
its terminals simultaneously, which is the property that we desire of 
variable elements. We do not wish to force any functional relation 
between the voltage across and the current through variable ele- 
ments. We wish only to observe constraints that may be imposed on 
the v — i relations by the VFE's and CFE's. The norator allows the 
network the extra degree of freedom taken away by the introduction 
of nullators. 

V. FORMULATION OF NETWORK EQUATIONS 

Since the introduction of nullators and norators into a network will 
generally introduce singularities into the corresponding equations, we 
call the approach we are considering the method of singular im- 
bedding. It has been demonstrated that the design problem can be 
reduced to the appropriate imbedding of nullators, norators, and in- 
dependent voltage and current sources. Let us now examine the effect 
of such imbedding on the network equilibrium equations. Since a 
nodal admittance formulation is used, it is important to determine 
the effect of nullators and norators on the admittance matrix. 

Independent voltage sources may be conveniently incorporated into 
an admittance formulation. If a series impedance exists with the 
voltage source, application of Norton's Theorem is sufficient. If no 
series impedance exists, the introduction of positive and negative im- 
pedances is necessary in transforming the voltage source to an inde- 
pendent current source (see Fig. 8) . 

The effect of nullators and norators upon the admittance matrix 
of a network has been considered by A. C. Davies. 5 Let us write the 
nodal equations for the network with all nullators and norators re- 
moved. The equations are of the form 

[Y„]V] = I,] (16) 



Fig. 7 — Schematic representation of norator. 
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Fig. 8 — Equivalent circuit for ideal voltage source. 

where 

[Y„] is the admittance matrix of the network with nullators and 

norators removed 
V] is the vector of node voltages with respect to ground 
I a ] is the vector of currents injected into each node. 

Suppose now that a nullator is connected between nodes i and ;'. 
Since the nullator passes only zero current, the current law equations 
at those nodes are not affected. However, since there is zero voltage 
across the nullator, V t and Vj are now constrained to be equal. Call 
this new value V {j . Clearly, one degree of freedom has been removed 
from the network response. In addition to the matrix equation (16), 
one equation of the form 

7, = 7, (17) 

is added for each nullator imbedded in the network. Thus, if k nulla- 
tors are imbedded, k additional constraint equations are added. 

Two viewpoints can be taken here. First, the original set of equa- 
tions, equation (16), has been appended by a set of the form 

[B]V] = 0] (18) 

where 

V] is the n-vector of node voltages. 

[B] is a k X n matrix of —1, 0, 1 entries expressing the set of con- 
straints of equation (17) for the k nullators. 

The final set of equations becomes 



ft] 



V] = 



(19) 



A- second approach to the problem was suggested by Davies. 6 In 
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the nodal equations below 
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(20) 



the addition of a nullator between nodes i and j makes V t = V, = Vy. 
The ith and ;th column of the Y matrix are both multiplied by Vy, 
thus they may be added and the equations written as 



Vu ••• (Vu + Vu) •" tfi. 

?y 21 (?/2, + Jhi) \hn 



V x 
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(21) 



The addition of k independent nullators (no nullator loops) causes k 
additions of columns of Y and reduces the dimension of V] by k. We 
denote the reduced set of equations by 

[Y' ] nxln - k) y'] ( „-„xi = I'mUi. (22) 

In either interpretation, we observe that the resulting set of equa- 
tions is no longer square. In the first interpretation, we are increasing 
the dimensionality of the vector space that the column vectors must 
span, without adding new basis vectors to span that space. In general, 
the equations will be inconsistent. In the second interpretation, we are 
keeping the dimension of the space fixed, but reducing the number of 
vectors available to form a basis and the space may no longer be 
spanned. Again inconsistencies will generally arise. In either interpre- 
tation, the inconsistencies are to be expected since nullators (VFE's 
or CFE's) have been introduced to constrain network variables. 

Let us now examine the way in which variable elements (additional 
degrees of freedom) remove these inconsistencies. Again two points of 
view may be taken. One provides us with new basis vectors to span 
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the space of possible injected current vectors I 5 ], the second reduces 
the dimensionality of the space of l s ] in order that the existing number 
of basis vectors might again span the space. 

Section III gives the essence of the first interpretation with the 
important result, equation (13). Observe that imbedding variable 
elements in a network provides an additional set of column vectors, 
namely, those of [ — C], that may be used as basis vectors in spanning 
the space of possible I a ]. Thus, if one has complete freedom in selecting 
variable elements, a set of column vectors, the columns of [ — C] can 
always be found to assure that the space of all possible I s ] will be 
spanned, regardless of how the nullators reduce the space of the column 
vector of the Y matrix. This concept, which involves growing new 
elements to satisfy imposed constraints, will be the subject of future 
study. 

A second approach in handling the freedom introduced by variable 
elements is to replace each variable element by a norator, as suggested 
in Section III. The method of Davies may then be employed to analyze 
the network containing norators. 6 Again assume that the admittance 
matrix Y„ of the network without nullators and norators is available. 
Thus, 

[Y„]V] = IJ. (23) 

Now suppose that a norator is connected between nodes h and k, and 
that the reference direction for the arbitrary norator current I is 
from h to k. The current-law equations for nodes h and k will be of 
the form 

I Sh ~ L - Z Y hiVi (24) 

l Sk + /. = L Y kiVi . (25) 

Since I is arbitrary, and is not needed to solve for the node volt- 
ages, adding the two equations gives 

/« + /« = E (Y hi + Y ki )v<. (26) 

This corresponds to the addition of rows h and k of the nodal equa- 
tions of the network without norators. Thus for a network containing 
n nodes and r norators, only n — 1 — r independent equations can be 
written. 

Observe in Fig. 9 that the effect of connecting the norator between 
nodes h and k is to replace the nodal equations for nodes h and k 
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Fig. 9 — Effect of connecting norator between two nodes. 



by a single current law equation for the ambit (broken line) sur- 
rounding both nodes h and k. Thus any functional relation between 
the current and voltage of branch ; is removed, as is desired for a 
variable element. 

To summarize thus far, the following manipulations may be per- 
formed on the network current law equations to deal with VFE's, CFE's 
and variable elements. To include network constraints, first imbed the 
CFE's and VFE's. Write the Y matrix with initiators removed. Then 
reduce the matrix by adding appropriate columns. This may be stated 
compactly by a matrix transformation as 5 



Is] = MUJV] 



(27) 



where [U c ] is a matrix obtained from the unit matrix by adding columns 
corresponding to nodes between which nullators are connected. Since 
the transformation [U e ] is singular, not all components of V] are deter- 
mined. The undetermined ones are found from the relation 



[B]V] = 0]. 



(28) 



To include variable elements, either 

(*) Augment the Y matrix of the fixed portion of the network with 
the node cutset matrix of the graph of the variable elements to get 
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= la] 



(29) 



or 



(m) Add the current law equation corresponding to nodes to which 
a nullator is connected. This is compactly stated by a matrix transforma- 
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tion as 



[U r ]I,] = [U r ][Y ]V] 



(30) 



where [U r ] is a matrix obtained from the unit matrix by adding rows 
corresponding to nodes between which norators are connected. The 
vector of currents through variable resistors is then formed by equation 
(29). 

VI. SOLVING THE NETWORK EQUATIONS 

We now wish to solve the set of equations after imbedding CFE's, 
VFE's, and variable elements. We assume equation (29) to be our 
starting point. A similar formulation may be made using equation 
(30) as the starting point. CFE's and VFE's are imbedded, variable 
elements are specified, and nullators are removed to generate the set 
of equations 
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Addition of nullators to the network adds the set of equations 

BV] = 0] 

and, from equation (27), the corresponding transformation [U e 
the admittance matrix. Thus the final set of equations becomes 

-c ! Y f U/ 
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(33) 



As seen in the previous section, the transformation [U e ] (which adds 
columns of Y f ) is consistent with the set of equations [B]V] = 0. Thus 
the second matrix equation in equation (33) will always have a solu- 
tion, provided the first one does. It remains only to solve 



-C ! Y f 



r fU c ]^] = I 5 ]. (34) 

in order to determine the proper element values. Let 

v] - * 

» J(r+n-*)Xl 

By using the Gauss- Jordan method one can bring these equations 
into the form 



L ! Oj X J I S2 J 



(35) 
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where 

Xi|X 2 ] is a vector of node voltages and currents through variable 

resistors, 
[U] is the unit matrix, 

[Q], I sl ], I S2 ] are the resulting submatrices after transformation. 
If I S2 ] = (the equations are consistent), the first equation can be 
solved for X x ] in terms of X 2 ]. 



X, = I 51 ] 



[Q]X S 



(36) 



The case I 52 ] 5^ implies that there are no values of variable elements 
consistent with the imposed constraints. For I S2 ] = 0], equation (36) 
generates all solutions to the problem. Some network variables X 2 ] 
can be chosen arbitrarily and the remaining variables X x ] determined. 
At each setting of X 2 ] the variable elements can be determined since 
all node voltages and currents through variable elements are known. 



Thus 



Zi = 






for 



i = 1 , r 



(37) 



where il and i2 are connection nodes of the ith. variable element. By 
allowing the free variables X 2 ] to take on a continuum of values, all 
solutions to the problem are determined directly. 

Returning to the example already discussed (Fig. b) , let us apply 
the method of singular imbedding. The circuit is redrawn in Fig. 10 
with the introduction of a VFE to constrain the voltage between nodes 
1 and 2 to be 0.1* volt. With the nullator removed, a set of nodal 
equations is written in the form of equation (13) 
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(38a) 



The introduction of a nullator between nodes 1 and 3 results in the 
addition of the corresponding columns and the equality 7i = 7 3 = 



* Since the nullator passes zero current, the series battery in the VFE model 
may have a nonzero resistance and still maintain the proper terminal voltage. 
Thus the introduction of positive and negative resistances are unnecessary here. 
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Fig. 10 — Network after singular imbedding. 



V 13 . Thus, 
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-1 
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= -0.1 
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v i3 
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(38b) 



With R = 1 ohm for ease of visualization, elementary row operations 
yield 



Thus, 



(39) 



(40) 
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o ! l" 


J, 


0.9" 


o i o ! i 
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= 0.9 


o o l ! -i_ 
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O.lj 


/," 


0.9" 




l" 




h 


= 0.9 


- v 2 


1 




v l3 . 


0.1 




-1 





It is clear that V 2 can take on arbitrary values while maintaining the 
constraints. We will demonstrate this for two particular values of 
y 2 .For7 2 = 
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111 


0.9] 


h 


_ 0.9 


v 13 


0.1 


v 2 _ 


0_ 



ffl = 



V, - v 2 



I, 



V, 



R 2 = y - °- 

M 

It is easily verified that these values, when substituting into the 
circuit of Fig. 1, result in V = 7, - V a = 0.1 volt. 
Similarly for V 2 = 0.6 volt 



111 


0.3] 


h 


_ 0.3 


v 13 


0.7 


v 2 _ 


0.6_ 



and R x = %, R 2 = 2. 

Again it is easily verified that V = Vi — V s — 0.1 volt. With this 
simple example in mind, let us consider the solution of more com- 
plicated networks by computer. 



VII. COMPUTER SOLUTION 

A program has been written to solve the design problem for resis- 
tive networks. The program performs the following operations 

(t) Accepts input of circuit description in conversational mode. The 
circuit may contain resistors (both fixed and variable), VFE's CFE's 
batteries, independent current sources, and current controlled current 
sources. 

(it) Generates C, Y f , and I 5 matrices for the network. 

(iit) Reduces equations to triangular form by a Gaussian reduction 
which pivots around largest elements in array. 

(iv) Those variables not in the basis after gaussian elimination are 
passed to the right side and stepped through specified range. Resistance 
values are printed for each setting of the free variables. Each set of 
resistance values will satisfy the given constraints. 

Four examples demonstrate the flexibility of the method. Suppose 
in the circuit of Fig. 11 one wishes to choose Rt and R 2 to provide 
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Fig. 11 — Transistor design problem. 

a collector current of 5 mA. A CFE of value 0.005 is placed in series 
with the collector and the circuit of Fig. 12 is fed into the program as 
in Table I. After the program sets up the equations and performs the 
gaussian elimination, it prints, that the voltage at node 3 is free. It 
can be arbitrarily chosen to generate sets of solutions. 

This free voltage is then, at the user's request, stepped from 7 volts 
to 10 volts in 1 volt increments. Combinations of Ri and R% which 
provide a collector current of 5 mA are printed in Table I. To verify 
these results the program dcanal 7 was used to determine the transistor 
collector current for the fifth set of resistor values. As the table shows, 
the collector current is 5 mA. 

A second example involves simultaneously constraining I = 5 mA 
and VCE = 5 volts. As Fig. 13 shows, R x , R 2 , and R 3 are variable. 
The network with a VFE and CFE imbedded is shown in Fig. 14, and 
the results given in Table II. Verification of the first set of resistance 
values is given. Observe that I c = 5 mA and VCE = 5 volts. 



12V -= 




R, 
Fig. 12 — Network after transistor modelling and singular imbedding. 



=- 12V 
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Table I — Printout of the Session to Solve the Circuit of Figure 11 

TYPE NO. OF BRANCHES, NODES, CONTROLLED SOURCES.BATTERIES, CURRENT SOURCES 
A=7 5 1 3 

TYPE BRANCH RESISTANCES 
B=l . I. 1260. 200. 1.E3 I .E4 I. 

TYPE FOR EACH BRANCH: INITIAL NODE, FINAL NODE, BATTERY NO. 
C=l 5 I 431 152 523 121 321 142 

TYPE VALUES OF BATTERIES 
D=0. 12. -.7 

TYPE FOR EACH CONTROLLED SOURCE! BRANCH NO. AND CONTROLLING BRANCH NO. 
E=6 4 

TYPE VALUES OF BETAS 
F=75. 

OPTION COMMANDS: DESIGN R 

TYPE NO. VAR. RESISTANCES, NO. VOLTAGE CONSTRAINTS, AND NO CURRENT CONSTRAINTS 

1=2 I 

TYPE BRANCH NO. OF VARIABLE RESISTANCES 

J=l 2 

TYPE BRANCH CURRENTS BEING CONSTRAINED 

M:6 

TYPE VALUE OF EACH CURRENT BEING CONSTRAINED 

N:.005 

THE FOLLOWING NODE VOLTAGES ARE FREE 
3 

ENTER LOWER LIMIT AND INCREMENT FOR EACH FREE VARIABLE AND NO. OF SETTINGS 
0=7. I . 4 

THE FREE VARIABLE: 7. 
R< l)=l .1850722E+O3 
R( 2>=9.9900O05E+02 

THE FREE VARIABLE: 8. 
R( I)=I.I84I104E+O3 
RC 2> =7. 9900003 E+02 

THE FREE VARIABLE: 9. 
R( 1)=! .I831496E+03 
R( 2)=5.9900OOOOE+O2 

THE FREE VARIABLE: 10. 
R( l):|.l821898E+03 
R( 2) =3. 9899998 E+02 

DESIGN COMMAND=KEEP 

ENTER VALUES OF FREE VARIABLES FOR DESIRED SET 

= 10. 

OPTION COMMAND: TRAN ALL 

VCE IC 

TRANS 9 

I 4.9398502 4.9999999E-03 




SET I c = 5mA 
Mr E = 5V 



Fig. 13 — Transistor design problem. 



CIRCUIT DESIGN 



295 



T 



12V 



1260 



751 >IOk 



200 — ► |(2J 

(D -AA/Hl|H 



VFE 5 



o 



1000 




© 

Fig. 14 — Network after transistor modelling and singular imbedding. 

A third example involves the rather complex three transistor circuit 
illustrated in Fig. 15. The imbedding of VFE's and CFE's to con- 
strain collector emitter voltages to 5 volts, and collector currents to 
10 mA is shown. 

Table III illustrates the results of a computer solution to the problem 
by the method of singular imbedding. Observe that currents through 
variable resistors 10 and 14 can be arbitrarily chosen and sets of resistors 
R 10 through R 1S generated. Four such sets are presented in Table III. 
Observe the results of an analysis indicating one such set properly 
biases the network. Table IV presents the results of an optimization 
program, based on pattern search, 1 to bias the network, for which 
forty-eight exploratory moves and 105 pattern moves were required. 
Each exploratory move involves between eight and 16 circuit analyses. 
Each pattern move involves an average of four analyses. Thus, approxi- 
mately 1000 matrix inversions are required. Since each inversion involves 
(n 3 )/3 operations, the number of operations to generate a single bias 
network ^243,000. 

Singular imbedding increases the number of nodes from 9 to 15. 
However, only one matrix inversion is required to generate a solution. 
Thus the number of operations =n/Z = 1125. 

Singular imbedding increases the efficiency in finding a solution to 
this problem by a factor of approximately 200. What is even more 
important is the ease with which equivalent networks are generated. 
Each equivalent network is generated by a matrix multiplication of 
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Table II — Printout of the Session to Solve the Circuit of Figure 13 

TYPE NO. OF BRANCHES, NODES, CONTROLLED SOURCES,BATTERIES, CURRENT SOURCES 

11:7 5 I 3 

TYPE BRANCH RESISTANCES 

B=l. I. 1. 1260. 200. I.E4 i. 

TYPE FOR EACH BRANCH: INITIAL NODE, FINAL NODE, BATTERY NO. 

C:l 3 I 121 541 132 3 2 5 421 152 

TYPE VALUES OF BATTERIES 

D=0. 12. -.7 

TYPE FOR EACH CONTROLLED SOURCE: BRANCH NO. AND CONTROLLING BRANCH NO. 

E=6 5 

TYPE VALUES OF BETAS 

F:75. 

OPTION COMMANDSrDESIGN R 

TYPE NO. VAR. RESISTANCES, NO. VOLTAGE CONSTRAINTS, AND NO CURRENT CONSTRAINTS 

1:3 I 1 

TYPE BRANCH NO. OF VARIABLE RESISTANCES 

J:l 2 3 

FOR EACH VOLTAGE CONSTRAINT, TYPE PLUS AND MINUS NODES 

K=4 2 

TYPE VALUE OF EACH VOLTAGE CONSTRAINT 

L=5. 

TYPE BRANCH CURRENTS BEING CONSTRAINED 

M=6 

TYPE VALUE OF EACH CURRENT BEING CONSTRAINED 

N=.005 

THE FOLLOWING NODE VOLTAGES ARE FREE 
3 

ENTER LOWER LIMIT AND INCREMENT FOR EACH FREE VARIABLE AND NO. OF SETTINGS 
0=4. 2. 2 

THE FREE VARIABLE: 4. 
R( l):6.36OI035E+02 
R( 2)=6.4979250E+02 
R( 3>=7.4140005E+02 

THE FREE VARIABLE=6. 
R( 1)=1. 2760788 E+03 
R( 2):|.0450373E+03 
RC 3):3. 4140001 E+02 

DESIGN COMMAND=KEEP 

ENTER VALUES OF FREE VARIABLES FOR DESIRED SET 

=4. 

OPTION COMMANDS:TRAN ALL 

VCE IC 

TRANS « 

I 5.0000E+00 5.0001 E-03 



the vector of free variables, which is stepped through a specified 
range, and the matrix of vectors not taken into the basis after tri- 
angulation. In this case the matrix is 19 X 2 and the vector of free 
variables is 2 X 1. Each multiplication involves 2 X 19 = 38 opera- 
tions. This means that up to 14,000 equivalent networks can be gen- 
erated with the same number of operations needed to give one solution 
by optimization techniques. 

The value of singular imbedding is apparent here. Only one equa- 
tion need be solved, and from it, all solutions are generated. 

As a fourth example, a network was designed for a specified %x\ 
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Fig. 15 — Three transistor network with VFE's and CFE's imbedded for desired 

biasing. 

and 2 2 i simultaneously. The circuit is given in Fig. 16. Ri,R 2 , and R 3 
are to be selected to give z n = % and z 21 = V 3 . After proper imbed- 
ding of VFE's and CFE's the network of Fig. 17 results. Table V 
gives the results of a computer run to design the circuit. The third 
set, Ri = R 2 = Rs = 2 is shown to give the desired z-parameters 
through the Y — A transformation of Fig. 18. 



VIII. RESISTOR CONSTRAINTS 

In many design problems it is desirable to constrain the values that 
the parameters take to lie within certain limits. For example, in 
biasing a transistor network, although solutions in which some resis- 
tors are negative are mathematically correct, in practice such net- 
works are unacceptable. 

If the designer has a good feeling for the circuit he is working 
with, his choice of the free variables resulting from gaussian elimina- 
tion with maximum pivoting will usually yield resistors with posi- 
tive values. There are, however, instances involving multiple feed- 
back paths where intuition cannot always be relied upon. In these 
instances it is possible that the values given by the designer to the 
free variables yield negative resistances. Furthermore, it may be 
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Table III — Printout of the Session to Solve the Circuit of 

Figure 15 

THE FOLLOWING BRANCH CURRENTS ARE FREE 
10 

14 

ENTER LOWER LIMIT AND INCREMENT FOR EACH FREE VARIABLE AND NO. OF SETTINGS 
0s|.E-3 I.E-3 i.E-2 .2E-2 2 

THE FREE VARIABLES ARE 

1.000E-03 I.000E-02 

R(lQ)=4.2366667E+03 
R(ll)=!.!028088E+09 
R(l2)=3.7O0R484E+03 
R(13)=4.0579244E+03 
R(14)=4.1685477E+03 
R(l6)r3.3704454E+03 
R(17)r3.78l I070E+03 
R(18)=4.l9l7495E+03 

THE FREE VARIABLES ARE 

2.000E-03 1.000E-O2 

R(I0>=2.1l83333E+03 
R<l t)=1.1372716E+09 
RCI2)=3. 9956586 E+03 
R(l3)=3.7284580E+03 
R(14):4.i685477E+03 
R(l6)r3.37O4454E+03 
R<17>=3.781 I070E+03 
R(IB):4. 1917495 E+03 

THE FREE VARIABLES ARE 

1.000E-O3 I.200E-02 

R<10>=4.2366667E+03 
R(ll):2.l 183293 E+03 
RU2)=3.7008484E+03 
R(13)=4.9290360E+03 
R<14)=3.4737897E+03 
R(l6):3.3704454E+03 
RC I 7) =3.781 1070 E+03 
R(I8>=4. 1917495 E+03 

THE FREE VARIABLES ARE 

2.0O0E-O3 1.200E-02 

RC 10> =2. I 183333E+03 
R(l l):2.H83293E+03 
R(12):3.9956586E+03 
R(l3)=4.45l26l5E+03 
R(l4)=3.4737897E+03 
RC16)=3.3704454E+03 
R(17)=3.781l070E+03 
R(18)=4.l917495E+03 

DESIGN COMMANDiKEEP 

ENTER VALUES OF FREE VARIABLES FOR DESIRED SET 

:2.E-3 I.E-2 

OPTION COMMANDSrTRAN ALL 

VCE IC 

TRANS # 

1 5.0000E+00 1.O000E-02 

2 5.O000E+00 1.0000E-02 

3 5.0000E+00 9.9999E-03 

difficult to explore the space of the free variables looking for regions 
where all the resistors are positive. 

One possibility for finding positive resistor regions is to use an 
optimization technique in which, considering the free variables as 
adjustable parameters, the sum of the absolute magnitudes of the 
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negative resistors is reduced to a minimum. If there exist solutions 
with all resistors positive, the minimum (zero) hopefully would be 
found automatically by the optimization routine. 

This optimization is more efficient than solving the problem by 
exploring a space in which all the variable resistors are parameters 
to be adjusted. 1 

Although the method given has been tried with success, a superior 
method having several advantages over the one proposed is explained 
in Section IX. The method avoids some of the most important problems 
associated with nonlinear programming. 

Some of these problems are : 

(i) The routine may get trapped in local minima. 

(u) Depending on the shapes of the surfaces involved and on the 
methods used the convergence towards the minimum may be very slow. 

(Hi) If the optimization is with constraints the nonlinear constraints 
are usually difficult to handle. 

If it were possible to reduce the problem to a linear programming 
problem, the following would have been gained: 

(i) If the problem has a finite minimum it will be achieved in a 



Table IV — Printout of Optimization Program 

INITIAL BRANCH RESISTANCES 

R(IO):0.5O00E+04 
R(l l)=O.50OOE+04 
R(l2)=O.3000E+04 
R(l3)=0.3000E+04 
RCI4>=0.3000E+04 
R<l6)=0.3000E+04 
RU7)=0.3000E+04 
RU8>=0.3000E+04 

EXPLORATORY MOVES 48 

PATTERN MOVES 105 

EINAL BRANCH RESISTANCES 

R(IO)=O.3730E+04 
R(l l)=0.4340E+04 
R( 12) =0.3732 Et-04 
RCI3):0.4377E+04 
R(l4)=0.3795E+04 
R( 16) =0.3372 E+04 
R(l7)=0.3783E+04 
R(l8)=0.4197E+04 

TRANSISTOR OPERATING POINTS 

TRANS # VCE IC 

1 5.000E+00 1.000E-02 

2 5.000E+00 1.000E-02 

3 5.000E+00 I.000E-02 
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Fig. 16 — Z-parameter design problem. 

finite number of steps. No local minima which are not also global 
minima exist. 

(ii) Algorithms exist which converge to the minimum efficiently. 

(Hi) The linear constraints generally complicate the problem only 
moderately. 

In Section IX the problem of biasing transistor networks is reduced 
to a linear programming problem. 



IX. APPLIED LINEAR PROGRAMMING 

Let us start by assuming a network in which the designer knows 
the correct signs of the node voltages with respect to the datum and 
the direction of the currents in the variable resistors. Generally the 
former is an easy task since it only involves knowing the nodes with 
the lowest potential. If this node is chosen as the datum, all the node 
voltages will be positive. Knowing the correct direction of the cur- 
rent through the variable resistors requires a better understanding of 
the circuit operation. Furthermore, there may be solutions in which 
the current through some resistors may flow in either direction. For 
this reason this requirement will eventually be relaxed. 

Linear programming requires the right side vector of equation (33) 

R 3 
V\A/ 



2/3 



2/3 



VFE(2/3) 



•2/3 



VFE(l/3) 



Fig. 17 — Network after singular imbedding 
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Table V — Printout of the Session to Solve the Circuit of Figure 16 

TYPE NO. OF BRANCHES.NODES, CONTROLLED SOURCES, BATTERIES, CURRENT SOURCES 

A:6 4 1 I 

TYPE BRANCH RESISTANCES 

B = l. I. 1. .66666667 .56666667 .66666667 

TYPE FOR EACH BRANCH: INITIAL NODE, FINAL NODE, BATTERY NO. 

C = l 2 I 141 241 231 341 311 

TYPE VALUES OF BATTERIES 

D=0. 12. -.7 

TYPE FOR EACH CONTROLLED SOURCEjBRAMCH NO. AND CONTROLLING BRANCH NO. 

E=6 5 

TYPE VALUES OF BETAS 

F:75. 

TYPE FOR EACH INDEPENDENT SOURCES INITIAL NODE AND FINAL NODE 

G:l 2 

TYPE VALUE OF EACH INDEPENDENT CURRENT SOURCE 

Hsl. 

OPTION COMMANDS=DESIGN R 

TYPE HO. VAR. RESISTANCES, NO. VOLTAGE CONSTRAINTS, AND NO CURRENT CONSTRAINTS 

Ir3 2 
TYPE BRANCH NO. OF VARIABLE RESISTANCES 
J = l 2 3 

FOR EACH VOLTAGE CONSTRAINT TYPE PLUS AND MINUS NODES 
X=2 1 4 1 

TYPE VALUE OF EACH VOLTAGE CONSTRAINT 
L:. 66666667 .33333333 

THE FOLLOWING BRANCH CURRENTS ARE FREE 

ENTER LOWER LIMIT AND INCREMENT FOR EACH FREE VARIABLE AND NO. OF SETTINGS 
Os-1. .33333333 3 

THE FREE VARIABLE = -1.0000 
R( l)= 0.66666670E+00 
R( 2>=-0.66666667E+00 
R( 3)=-0.66666680E+00 

THE FREE VARIABLE = -0.6667E+00 
IK I): 0.1 0000000 E+0 I 
R( 2>r-0.20000000E+0l 
R( 3>s-0.l9999999E+0l 

THE FRE VARIABLE = -0.3333E+OO 
R( 1): 0.I9999997E+01 
R( 2): 0.20000000E+0I 
R( 3): 0.2000001 3E+0 I 
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Fig. 18 — Verification of computer solution. 
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to have positive entries. This may be achieved by multiplying by —1 
all those rows in equation (33) which have a negative entry in the 
right side vector and thus obtain the set of equations 



(41) 



_0 


H 2 " 
H 3 _ 


I] 
V 


f] 
0_ 



where 



is obtained from 



H = 



_0 


H 2 ] 
H 3 _ 



r -C 


Y f U c " 


f 

and - 



_ 


B 



from — 



by possibly multiplying some rows by —1. 

To force all the branch voltages to be positive let us add the constraint 

-[C]V] ^ 0] (42) 

where [C] is the matrix appearing in equation (12). 

Equation (41) and inequality (42) together with the condition 

£] ^ 0] (43) 

can be looked upon as a linear programming problem in which it is 
desired to find the value of a positive vector satisfying a set of linear 
equalities and inequalities and which minimizes the linear function 
where 



D = [0, 0, • • • , 0]. 



(44) 



Since the minimization of the constant zero is of no interest, all that 
is required is to obtain the feasible solutions of the linear programming 
problem. 6 

Once the feasible solutions are obtained, the fact that the solution 
satisfies equation (41) guarantees that the circuit is properly biased 
while the positivity condition on the vectors I|V and — [C]V] guarantee 
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that all the variable resistors are positive, since both the currents 
and voltages across them are positive. 

To obtain the feasible solutions phase I of the two phase simplex 
method may be used. 6 

Phase I of the simplex method finds the basic positive solutions of the 
system of equations 



H, 


H 2 





I" 

V 


f] 
= 
0_ 





H 3 





Lo 


-C 


-u 



(45) 



where the vector w] (which is constrained to be positive) is a slack 
vector and U is a unit matrix. 

By denoting with A the matrix on the left of equation (45), with x 
the column on the left, and with b the column on the right side, equa- 
tion (45) may be written 

Ax = b. (46) 

Let the dimensions be: A, m X n;x, n X 1; b, m X 1. Let A and [A | b] 
have rank r. This implies equation (46) is compatible. (The case in 
which this is not true is of no interest since in such case no solution — 
whether positive or not — exists.) 

Phase I of the simplex method finds positive solutions of equation 
(46) for r of the variables x i} i — 1, 2, ■ • • , r setting the rest of the 
x } , j = r + 1, • • • , n to zero." Each one of this set is a basic feasible 
solution. There may be several such sets for a given problem. The 
totality of nonnegative solutions of equation (46) is the convex hull 
of the basic solutions. By extending the simplex algorithm so that 
once a basic feasible solution is found the other basic feasible solu- 
tions are also searched for, it is possible to obtain all basic feasible 
solutions. 



Suppose x l , x 2 , 
x satisfying 



x p are basic feasible solutions. Then any vector 



x = X.x 1 + X 2 .r 2 + • • • + A/..r' 
X 1 , X 2 , • ■ • , \ P ^ 

x, + x 2 + • • ■ + X P = 1 
is also a feasible solution. 



with 
and 



(47) 



* In case no nonnegative solutions to equation (46) exist, the simplex algorithm 
is able to detect it. 
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If x 1 , x 2 , • • , x p is the set of all basic feasible solutions, then all the 
solutions of equation (47) constitute the complete set of feasible solu- 
tions. 



X. RESISTORS WITH UPPER AND LOWER BOUNDS 

In the previous discussion the appearance of nonnegative resistors 
was precluded by adding inequality (42) . Often it is desirable to im- 
pose lower and upper bounds for the resistors because the technology 
used to realize them requires it. For example, if tantalum thin film 
resistors are used it is desirable to restrict them to lie between 10 and 
10 fl ohms. 

Let the fcth variable resistor be connected from node i to node ;. The 
value of R k is given by 



R* = 



7,- - 7, 



If it is desired to have this resistor lie within 10 and 10" ohms the 
following conditions are imposed 



^ 10, hi* 



7, 


— 


V, 




h 




Vi 


— 


v. 



h 



^ 10 5 , I k 9*0 



(48) 



which may be rewritten (recall I k is nonnegative) 
7, - 7, - 107* ^ 

7, - 7, - 10% ^ 

If instead of equation (42) inequalities similar to equation (48) 
are written for all variable resistors, the resulting circuits will have 
all variable resistors within specified upper and lower bounds (except 
for the possibility I k = 0, which implies an open circuit, in which 
case the resistor disappears altogether) . 

The problem of biasing of transistor networks with positive resis- 
tors is equivalent to solving 



_o 


H 2 
H 3 _ 







IM^]*Sl 



(49) 
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[D 2 



-*-i\*i 



.i»a 



where D! and D 2 are diagonal matrices whose diagonal elements contain 
the minima and maxima for the variable resistors. By adding positive 
slack vectors Wi and w 2 , equation (49) is equivalent to 
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_D 2 


-C 





uj 



(50) 



where U is a unit matrix and the vector on the left is restricted to be 
nonnegative. 



XI. RELAXING SIGN CONDITIONS 

So far it has been assumed that the direction of the current flow in 
variable resistors is known beforehand. This condition may not hold 
for some cases and hence it is desirable to relax it. 

When a variable in a linear programming problem is not required 
to be positive it is customary to write it as the difference of two posi- 
tive quantities. Thus if I k and Vj - V k = V are not required to be 
positive one may write 

h = /*• - h- 

Vi = v l . - v r . 

where I k . , I„. , V v , V v . ^ 0. 

A current I k of variable sign may be restricted to have a magnitude no 
less than I ok ^ by imposing the pair of conditions* 



h ^ L 



or 



■/* ^ /. 



(51) 



Likewise a branch voltage V t of variable sign across a resistor may be 

* The constraint set on the currents is not convex, therefore it is necessary to 
solve the problem twice, once with each inequality, and take the union of the 
two solutions. If n variable resistors may have currents flowing in either direc- 
tion, the solution will be the union of' the solutions of 2" problems in which 
all the combinations of the inequalities are used. 
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restricted to have a magnitude no greater than V ok ^ by imposing the 
pair of conditions 

V t ^ V ok and -V l ^ V ok . (52) 

If Vi is the voltage across the kth. resistor and I k its current, then in- 
equalities (51) and (52) insure that the magnitude of the kth resistor 
satisfies 

V < * 
\Rk\ U Y"- (53) 

The resistor R t may be negative or positive. However, if each variable 
resistor is made of two resistors in series one of value V„ k /I ok and the 
second to be determined by the computer subject to equation (53), the 
series combination of the two resistors will never be negative. This 
constitutes a technique for guaranteeing positive variable resistors 
without previous knowledge of the directions of current flows.** 

The method described can also handle circuits with variable resistors 
whose values he within upper and lower limits. If R kml „ and R kmai are 
the minimum and maximum values allowed for the /cth variable resistor, 
the fixed series resistor should be 

R k , = R kmin + VJI ok (54) 

with V Bk /I ok chosen such that 

V ok /I ok = R km „ - R kmin . (55) 

The value of R km , „ may be zero. Thus, a resistor may disappear as a short 
circuit. If instead of bounding the value of a resistance from above, the 
value of an admittance is bounded, a dual method may be used to 
guarantee positive resistors. 

Instead of equations (51) and (52) the following restrictions are 
imposed 

7, ^ I ok and -I k % I ok , (56) 

V, ^ V ok or -Vi ^ V ok A (57) 



* Roth V.h and /.* are variables in the linear program which will be determined 
by the simplex algorithm. The ratio is constrained by a linear inequality hk \ Rk | 
—Vok < 0, where \ R*\ is given. 

** Another approach is to reverse the reference direction of the current and 
voltage drop across each variable resistor and apply the methods of the previous 
section. If n variable resistors may have currents flowing in either direction 
it is necessary to consider 2" possibilities. 

t See footnote to equation 51. 



CIRCUIT DESIGN 307 

These guarantee that 

|0»l ^ 4r (5 ® 

I ok 

where G k = \/R k . G k may be positive or negative if each variable 
resistor is made of two resistors in parallel, one of admittance I„ k /V ok 
and the second to be determined by the computer subject to equation 
(58). However, the parallel combination of the resistors will never be 
negative. 

The dual method can also handle circuits with variable resistors whose 
admittance lies within upper and lower limits G km ^ and G kmin . The value 
of G kmln may be zero. Thus a resistor may disappear as an open circuit. 

XII. CHOOSING TOPOLOGY BY COMPUTER 

As already pointed out, Phase I of the simplex method obtains the 
basic feasible solutions of a set of linear equations. The set of equa- 
tions may come from a set of equalities and inequalities to which 
slack variables have been added. Usually the number of variables 
(including slack variables) is greater than the number of equations 
and the system is redundant. If r is the rank of the system and n is 
the number of variables (including slack variables) , at least n-r 
variables are set to zero in obtaining a basic feasible solution. Some of 
the variables set to zero may be node voltages or variable resistor 
currents. If a node voltage is set to zero, the corresponding node is 
grounded. If a variable resistor current is set to zero, the corre- 
sponding resistor disappears as an open circuit. If a slack variable 
is set to zero, the inequality constraints are met with equalities. 

For example, for equation (50) if the fcth entry of w, is zero, the 
fcth resistor acquires its minimum allowed value. 

One way of viewing equation (50) is to consider the columns of 
the matrix on the left as elements of a vector space and the entries 
of the column multiplying the matrix as those positive coefficients 
which synthesize the column on the right in the form of a linear com- 
bination of the columns of the matrix. A final tableau of Phase I of 
the simplex method will contain a number of independent unit col- 
umns (with all entries zero except one) equal to the rank of the mat- 
rix on the left side of equation (50). The unit columns are obtained 
by the special gaussian reduction provided by the simplex algorithm. 
Each column corresponds to a variable in the column multiplying the 
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matrix of equation (50). Those variables whose corresponding col- 
umns are not unit columns are set to zero. 

If a set of columns corresponding to the currents through a set of 
variable resistors are linearly dependent, one or more of the currents 
will be set to zero. This implies the disappearance of a resistor as an 
open circuit. The choice of which resistors disappear is automatically 
determined with the aid of the simplex algorithm, so that the non- 
zero currents acquire positive values (if such a choice exists). If two 
columns of the matrix of equation (50), corresponding to currents 
through variable resistors, are linearly dependent it means that Kirch- 
hoff's voltage and current law may be satisfied with one of the cur- 
rents zero, making one of the resistors unnecessary. 

The above argument provides a method for letting a computer pro- 
gram choose the topology and resistor values of a dc network in 
which certain voltages and currents are imposed by CFE's and VFE's. 
One connects an excess of resistors between different nodes (includ- 
ing additional internal nodes if desired) . By using a linear program- 
ming formulation some node voltages and variable resistor currents 
are set to zero by the computer program, thus determining a set of 
"linearly independent positive resistors" that satisfy all the circuit 
equations. 

XIII. EXAMPLES 

Consider the circuit of Fig. 19(a). The equivalent circuit is shown 
in Fig. 19(b) with a VFE and CFE in place. As indicated on Fig. 
19(b) it is desired to impose on the transistor a collector current of 
5 mA and a collector-emitter voltage of 5 volts. The resistors marked 
Ri, R 2 and #3 are variable. 

The nodal equations for the circuit after the effect of the nullators 
introduced by the VFE's and CFE's are taken into consideration 
are, in matrix form 



r -l. 


1. 


0. 


0.005 




-0.005 


0. 


-1 X 10~ 10 


0. 


0. 


1. 


-0.38 




0.3811 


-0.0011 


0. 


0. 


0. 


0. 


0. 




0. 


0.00333 


0. 


1. 


0. 


0. 


-1. X 


1Q -io 


0. 


0. 


1. 


0. 


0. 


0. 


0. 




-0.001 


0.001 


0. 


0. 


0. 


0. 


0.375 




-0.3751 


0.0001 


0. 
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(a) CIRCUIT 



(b) MODEL 



Fig. 19 — Circuit biased with constrained singular imbedding; (a) circuit, (b) 
model. 



L 




0.0035 


h 




-0.27 


h 




0.035 


v 2 


= 


12. 


v 3 




0.005 






0.2675 



(59a) 



By multiplying the second row by —1, the entry —0.27 in the right 
side vector is made positive. (As indicated above, linear program- 
ming assumes the right side vector is nonnegative) . Notice that since 
the matrix in equation (59a) is 6 X 7, we therefore generally expect 
a one parameter infinity of solutions. If the system of equations were 
solved using the simplex method (with arbitrary cost coefficients) , 
solutions in which all the variables acquire non-negative values may 
be obtained. Resistors R 2 and R s , which are grounded, will auto- 
matically be positive. However the voltage differences across un- 
grounded resistors may turn out to be negative, yielding negative re- 
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sistanccs. To assure a non-negative voltage difference across R, the 
following additional constraint will be imposed 

V 5 - V 2 ^ 
which may also be written* 

V 2 - V 5 =g 0. (59b) 

There arc two basic feasible solutions to this problem: 



u 
I, 

V 2 
V 3 

v< 



5.788 

5.78794 

5.060073 X 10~ 3 
6.212 
5.5 
10.5 
6.212 



5.9966425 X 10" 5 

0. 

5.060073 X 10" 3 

6.212 

5.5 
10.5 
11.99994 



The first basic solution yields the following set of resistors 



Ri = 



V, - v 2 



6.212 - 6.212 



u 



6.212 



a * ' U 5.78794 



/? 3 - 



W 



5.788 
= 1 .073266 ohms 
5.5 



= ohms 



35 = 1086.941 ohms. 



(60) 



(61) 



U 5.060054 X 10 
/?, is a short circuit. 

The second basic solution yields the set 

R x = 97008.514, R 2 = », 7^ 3 = 1086.941. 

R 2 is an open circuit. Notice also that R 3 is the same for both solu- 
tions. This is expected since the voltage of node 3 is virtually fixed 
by the requirements. 

The totality of the solutions with non-negative voltage differences 
across the variable resistors may be written, according to equation (47) 

x = Xx 1 + (1 - \)x 2 

where x l and x 2 are the basic feasible solutions of equation (60), and 
< X ^ 1. 



* When the right side of an inequality is zero, it is preferable to write it as 
a ^ inequality because the corresponding slack variable may be used as an 
artificial variable with savings on the size of the matrix to be manipulated. 
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Choosing X = \ yields 



x = 



5.5 
10.5 

. 9.10597 
which yields the set of resistors 

Rx = 0.99790G, R 2 = 2.1-16532, 



2.89402 
2.89397 

5.060073 X 10" 3 
6.212 



H, 



1086941. 



A continuous set of equivalent circuits, which achieve the require- 
ments exactly and which have positive resistances, is obtained by 
varying A between to 1. 

Suppose now that further considerations require that R x lie be- 
tween 1000 and 2000 ohms. By replacing (59b) by 



V 6 - v a 



> 1000 



and 



V, - V, 



< 2000 



which may be written 

1000 /, - F 5 + V 2 ^ 
-2000 J, + V s - V 2 ^ 
the resistor R x is forced to remain between 1000 and 2000 ohms. 



(62) 



When the new problem is solved the basic feasible solutions are 

5.782218 X 10" 3 " 
5.722218 X 10" 3 
5.060073 X 10" 3 



x = 



6.212 

5.5 
10.5 
11.99422 



2.892554 X 10" 3 
2.832554 X 10~ 3 
5.060073 X 10" 3 
x 2 = 6.212 
5.5 

10.5 

11.99711 
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The basic feasible solutions yield the following sets of resistors 
x x :R { = 1000., R 2 = 1085.593, R 3 = 1086.941 
x 2 :#, = 2000., R 2 = 2193.074, R 3 = 1086.941. 

Notice that R t acquired its allowable extreme values in each basic 
feasible solution. 

Other sets of resistances may be obtained by convex combinations 
of the two basic feasible solutions. 

As an example in which the topology of a circuit is determined by 
the computer, consider the circuit of Fig. 16 in which B lf R 2 , and R s 
are to be selected to give z lt = % and z 21 = %. The example was 
previously solved without linear programming techniques. Several 
solutions appear in Table V. By maximizing the negatives of the cur- 
rents in the resistors, those currents which may be set to zero by tak- 
ing them out of the basis for a basic feasible solution will be con- 
verted into open circuits. After the effect of the nullators introduced 
by the VFE's is accounted for, the matrix corresponding to the 
circuit of Fig. 16 is 5 X 6. We therefore expect a one parameter in- 
finity of solutions and two basic feasible solutions which are 






/l~ 


0.0 


0.28174743 


h 


5.9652404 X 10" 3 


4.5937138 X 10" 


u 


5.965238 X 10 " 3 


0.0 

, x = 


v 2 


44.096881 


43.602828 


v 3 


0.3333333 


0.33639577 


v 4 _ 


0.3333333 


0.3333333 



Notice that V 4 remains constant for both basic feasible solutions. This 
is expected since a VFE is connected from node 4 to node 1 (datum) . 
The resistances corresponding to the basic feasible solutions are 

a- 1 :R X = oo R 2 = 55.879273, R 3 = 7336.429 

x 2 : R, = 154.4759, R 2 = 72.5629, R 3 = «> . 

In both basic feasible solutions one of the resistances disappeared as 
an open circuit. This indicates that given Ri, R 5 and R Q with the 
values indicated in Fig. 16 the circuit is achievable with two topologies, 
each containing 5 resistors. 

Let us now make R 9 a variable resistor. The nodal matrix after the 
elimination of the nullators is now 5x7. Thus we expect a two 



CIRCUIT DESIGN 313 

parameter infinity of solutions and at least 3 basic feasible solutions. 
The following sets of resistors correspond to basic feasible solutions 

(i) R t = oo, R 2 = oo, R 3 = 1.33333, R e = 0.16666667 

(ii) R, = oo, R 2 = 0.333333, R 3 = 0.4444456, R a = co 

(iii) R, = 0.88888898, R 2 = 1.333333, R 3 = co , R 6 = oo 

(iv) R, = 1.3333336, R 2 = °o , R a = », # 6 = 0.66666651. 

These sets provide four different topologies with which given two 
of the resistors (72 4 and R$) a resistive network having z n = %, z 2 i 
= % may be realized. 

The example illustrates how using the methods of this paper can 
solve the problem of realizing portions of a resistive matrix with cer- 
tain elements prespecified. The prespecified elements need not be 
resistors but may also include controlled sources, gyrators, ideal trans- 
formers, and so on. 

The methods discussed have been implemented on a time-shared 

Table VI — Printout of the Session to Solve the Circuit of 

Figure 19 

TYPE NO. OF BRANCHES, NODES, CONTROLLED SOURCES, BATTERIES, CURRENT SOURCES 

A:7 5 I 3 

TYPE BRANCH RESISTANCES 

•iv I . I. I. 200. I.E4 300. 1. 

TYPE FOR EACH BRANCHj INITIAL NODE, FINAL NODE, BATTERY NO. 

C:5 2 I 211 311 232 431 143 153 

TYPE VALUES OF BATTERIES 

D:0. -.7 12. 

TYPE FOS EACH CONTROLLED SOURCE: BRANCH NO. AND CONTROLLING BRANCH NO. 

E=6 5 

TYPE VALUES OF BETAS 

F=75. 

OPTION C0ni1AHDS = DESIGN CKT _«..«.- 

TYPE NO VARIABLE RESISTANCES, NO. VOLTAGE CONSTRAINTS, AND NO CURRENT CONSTRAINTS 

1:3 1 I 
TYPE BRANCH NO. OF VARIABLE RESISTANCES 
J:l 2 3 

TYPE PLUS AND MINUS NODES FOR EACH VFE 
K=4 3 

TYPE VALUE OF EACH VFE 
L:5. 

TYPE BRANCH CURRENT FOR EACH CFE 
"1:5 

TYPE VALUE OF EACH CFE 
N:.O05 

TYPE COST COEFFICIENTS 
0=1. t. 1. 1. I. I. I. 
TYPE MINIMA OF VARIABLE RESISTANCES 
P:|0a0. 3. 0. 

TYPE MAXIMA FOR EACH VARIABLE RESISTANCE 
Q:200O. I.E8 1.E8 

R( 1): 9.9999995E+02 
R( 2): I .085593OE+03 
RC 3): 1. 0869498 E+-03 
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computer system. The program is conversational. A portion of a ses- 
sion in which a basic solution corresponding to the circuit of Fig. 19 
with Ei constrained betwen 1000 and 2000 ohms appears in Table VI. 

XIV. CONCLUSIONS 

The method of singular imbedding has been shown to be efficient for 
solving the following problem: Given a circuit with a prespecified 
topology, some of whose elements are prespecified, find the values of 
the unspecified elements which will yield desired node-pair voltages 
or branch currents. The unspecified element values may be restricted 
to lie within given upper and lower bounds. 

By letting the upper and lower bounds become infinite and zero, 
the problem of finding the topology for the circuit may be also solved. 

The method has been implemented on a time-shared computer, 
and several examples, including some practical transistor circuits, are 
given. 

The usual approaches to the problems of this paper have been itera- 
tive analysis-optimization schemes. Singular imbedding requires, for 
a three transistor amplifier, three orders of magnitude less computa- 
tion time. This makes the method appealing for time-shared applica- 
tions. 

Two new singular network elements, the voltage forcing element 
and the current forcing element, constrain node-pair voltages and 
branch currents without otherwise affecting the circuit. Elements of 
unspecified value are modeled by branches carrying unknown cur- 
rents. 

With the aid of these elements, the problem of design is reduced to 
one of analyzing a circuit containing unknown current sources and 
nullators. If there are more free elements than requirements, the solu- 
tion space may be a linear manifold. By allowing the free circuit 
variables to take on a set of discrete values, sets of exact solutions 
to the design problem may be generated economically. 

When the unspecified elements are required to lie within upper and 
lower bounds, the problem is one of analysis with linear inequality 
constraints. This may be solved efficiently using linear programming 
techniques. 

Among the practical problems solved by singular imbedding are 
biasing a direct coupled transistor amplifier, designing midband gain 
and driving point impedance, synthesizing networks for several given 
admittance parameters, and determining circuit topology. 

Areas being investigated include using singular imbedding in the 
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synthesis of resistance networks (the synthesis of a single column of 
a specified resistance matrix has been illustrated). Synthesis of an 
entire resistance matrix results from the intersection in resistance 
space of the solution spaces for each column of the matrix. Similarly, 
by considering the intersection of solutions spaces for both a small 
signal design and a biasing design, the method may be extended to 
designing transistor circuits for desired small signal design and bias 
points simultaneously. 

Although only fixed value CFE's and VFE's were used in this 
paper, CFE's and VFE's which may take any value within a given 
range may also be used. For example, a branch current may be forced 
to be greater than 1 mA and less than 10 mA. These elements are also 
useful in insuring that models for devices stay within their valid 
limits. For example, a transistor can be constrained to remain in the 
active region, for which the linear model used is valid. 

For simplicity, only the case of linear dc networks has been il- 
lustrated in this paper. However, the method has usefulness in ac 
design, combined ac and dc design, and non-linear design. These topics 
will be covered elsewhere. 
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